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Theoretical studies on the geometry, electronic structure and spin multiplicity of Sc, Ti and V 
doped Na n (n = 4, 5, 6) clusters have been carried out within a gradient corrected density functional 
approach. Two complementary approaches including all-electron calculations on free clusters, and 
supercell calculations using planewave pseudopotential and projector augmented wave formalisms 
^ , have been carried out. It is shown that spin magnetic moments of the transition metal atoms, the 

magnitude of host polarization, and the sign of the host polarization all change with the number 
of alkali atoms. In particular the transition metal atoms are shown to attain spin moments that 
are higher than their atomic values. The role of hybridization between the transition atom d-states 
and the alkali sp-states is highlighted to account for the evolutions in the spin moments and host 
polarization. 
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I. INTRODUCTION 



Extensive research over the past two decades has shown that atomic clusters consisting of two to a few hundred 
atoms constitute a new phase of matter with distinct properties different from those of individual atoms on the 
one hand and bulk matter on the other— Of all the properties, the magnetic properties in reduced sizes are 
probably the most fascinating 7 . Small clusters of itinerant ferromagnetic solids Fe, Co, and Ni are found to possess 
enhanced magnetic moments that can be around 30% higher than the bulk moments^. These clusters also exhibit super 
paramagnetic behavior as a consequence of reduced magnetic anisotropy as opposed to the ferromagnetic behavior 
observed in solids&2Jii. Small clusters of non-magnetic solids like Rh are found to be magnetic with appreciable 
moments per atom that disappear as the size approaches around 100 atom o n i 12 . Transition metal (TM) atoms doped 
in clusters of other metals also show very interesting magnetic properties. Recently, Li et al. — studied 3d TM doped 
i i ' Aug clusters both experimentally, and theoretically through density functional theory (DFT) calculations. In these 
clusters, the TM atoms were found to retain their atomic moments. Later, Torres et al. — reported DFT calculations 
of TM doped small cationic Au clusters in which the magnetic moment showed pronounced odd-even effects as a 
function of the cluster size, and resulted in values very sensitive to the geometrical environment. In all these cases, 
the moment per atom in the clusters are less than or equal to the possible spin magnetic moment in the individual 
atom. It has however been recently reported that atomic and even magnetic moments beyond the atomic value on 
transition metal atoms can be stabilized by depositing them on alkali metal surfaces and in the interior of alkali hosts. 
For example, Bergmann and Son g 15 i 16 have carried out anomalous Hall effect measurements on systems consisting of 
Vanadium impurities on Na and K thin films. Vanadium is non-magnetic in the solid phase. Yet, the measurements by 
j — ' Bergman and Song indicate a large moment in the range of 6.6 to 7 [Ib per atom. Sahu and Kleinmani 7 - have carried 
out first principles theoretical studies on Vanadium inside alkali hosts using a supercell approach. The theoretical 
studies suggest that the absolute value of the moment may be less than that predicted by Bergman and Song. Both 
studies, however, do confirm large spin moments on V sites indicating that alkali hosts can stabilize and enhance spin 
moments on transition metal sites. 

In this paper, we present an intriguing finding where the spin magnetic moment of transition metal atom and the 
polarizability of the host atoms can be controlled by changing the size of the alkali host. Our investigations cover 
Sc, Ti, and V atoms deposited on alkali templates consisting of four, five and six sodium atoms. The investigations 
employ first principles theoretical studies within a gradient corrected density functional to incorporate the exchange 
and correlation effects. We demonstrate that as the cluster size is changed, the three impurities exhibit differing 
behavior. While the spin moment of Sc and V atoms oscillate as the number of alkali atoms is changed, the moment 
of Ti atom monotonically increases as the size is increased. A detailed analysis of the atomic spin charges on Na 
atoms shows that while the host polarization enhances the spin moment in some cases, it opposes that in others. The 
ability to control the spin moment and the host polarization by changing the metal atom and the size of the alkali 
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cluster offers unique opportunities for tuning magnetic moment and coupling in such systems. 

II. COMPUTATIONAL METHODS 

All the theoretical studies were carried out within the gradient corrected DFT framework. We carried out all- 
electron calculations with atom-centered Gaussian type orbitals using the deMon2K codei 8 -. These calculations 
employed the PW86 generalized gradient approximation (GGA) functional^, and the DZVP basis sets optimized 
for gradient corrected exchange-correlation functional 20 . In order to avoid the calculation of four-center electron 
repulsion integrals, the variational fitting of the Coulomb potentia l 21 ' 22 was employed. The auxiliary density was 
expanded in primitive Hermite Gaussian functions using the A2 auxiliary function selj 2 ^. The exchange-correlation 
potential was calculated with the orbital density. Structures were optimized using the quasi-Newton method in internal 
redundant coordinates 2 ^ for different fixed spin multiplicities, and without any symmetry constraints. Supplementary 
calculations were performed in order to eliminate any uncertainty due to the basis set or the numerical procedure. 
We also carried out plane- wave pseudopotential calculations using a supercell approach. Here, an energy cutoff of 500 
eV was used, and the cluster was placed in a large cubic box of sides 20 A in order to reduce its interaction with 
its images. The potential between the ion cores and the valence electrons were expressed in terms of the Vanderbilt- 
type ultra-soft pseudopotential (USPP) 25 . The exchange-correlation effects were treated with the PW91— GGA 
functional. Brillouin zone integrations were carried out using only the T-point. Structures were relaxed using the 
conjugate gradient (CG) method for different fixed spin multiplicities, and without any symmetry constraints. We 
re-calculate total energies of all clusters using projector augmented wave (PAW) method in the same geometry as 
obtained from USPP calculations. These calculations were carried out using VASP cod o 27 ' 28 ' 29 . It turned out that the 
all-electron, USPP and PAW calculations all agreed on the lowest energy spin multiplicity for each cluster. We found 
some differences in the energy orderings and spin multiplicities of excited states between the plane- wave and Gaussian 
type methods; however these differences are small in most of the cases. In what follows we report the geometries and 
energy values based on all-electron deMon2K calculations, unless otherwise mentioned. 

III. RESULTS AND DISCUSSION 

A. Pure Na clusters 

In order to benchmark our calculations, we first calculated structures and binding energies (BE) of pure Na ra clusters 
in the size range n = 3, 4, . . . , 10 where previous calculations exist. Binding energy of a cluster is defined as the energy 
gain in assembling the cluster from the isolated atoms. For pure Na clusters this is given by the equation: 

BE(N& n ) = n£(Na) - £(Na„), (1) 

where E(Na„) represents the total energy of a Na„ cluster, and E(Na) is the total energy of an isolated Na atom. 
Note that according to this definition, BE is a positive number for a bound structure, and a larger BE implies a more 
stable structure. 

In both all-electron and USPP methods we were able to reproduce the previously reported structures of pure Na 
clusters 30 . In particular, we reproduced the C2„ structure for a pure Nas cluster, and for a Na6 cluster we found two 
close lying isomers with Cs„ and D^h symmetries. We found the 3-dimensional (3D) Cs„ structure to be marginally 
(~ 0.03 eV) lower in energy than the planar D 3 h structure. Earlier calculation s 30 ! 31 also found these two structures 
to be very close in energy with different methods giving different energy ordering. We also found that the BE per 
Na atom shows a peak at n = 8 relative to neighboring sizes. Enhanced stability of Nas with 8 valence electrons 
has been previously rationalized within a jellium model^ 2 . where the cluster corresponds to a magic species with filled 
electronic shell. 

Having confirmed that our calculations reproduce the known structures and energetics of pure Na clusters at small 
sizes correctly, we then studied TM doping in Na„ clusters for n = 4 — 6. At each size, we tried several initial 
geometries and relaxed them for different fixed spin states without any symmetry constraints. We now discuss the 
results of these calculations. Figure [1] presents the optimized geometries and Mullikcn atomic spin charges of the 
ground state and the next higher lying isomer in energy for the series TMNa„ (n = 4, 5, and 6 and TM = Sc, Ti, 
and V), and in Table [J we present the magnetic moment and relative energies of these clusters. 
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B. TMNa 4 clusters 

In the case of TMNa4 clusters, we found two low lying structures for ScNa4 which are shown in Figure [TJ We found 
the same structure for the ground state for ScNa4, TiNa4 and VNa4 clusters. This structure is obtained by replacing 
the Na atom with the highest coordination number in the C27, structure of a pure Na5 cluster by a TM atom. ScNa4 
and TiNa4 clusters retain the atomic moments of the Sc and Ti atoms, 1 and 2 fis respectively. What is intriguing 
and amazing is that the VNa4 cluster has an enhanced magnetic moment of 5 us, while that of a V atom in the 
ground state is only 3 (1b- The microscopic reason behind this enhanced magnetic moment is discussed later. 

A summary of the results on these clusters are presented in Table |U It is to be noted, however, that both for TiNa4 
and VNa4 clusters, the high and low spin-multiplicity states are very close in energy, while the two structural isomers 
of ScNa4 are also close in energy. 

C. TMNa 5 clusters 

For TMNas clusters, the low energy structures we found are: (i) a Cs v structure with a TM atom incorporated at 
the center of a Nas pentagon for the ground states of TiNas and VNas; (ii) a square bipyramidal structure having C4„ 
symmetry for the ground state of ScNa5 and the low lying isomer of VNas and, (iii) a structure with C s symmetry 
that is obtained by distorting the C4„ structure for the low lying isomers of ScNa5 and TiNas ■ These structures are 
shown in Figure [1] and the main results of our calculations for TMNa 5 clusters are presented in Table |U 

Although the ScNas cluster has no net moment in its lowest energy structure, we found a nearly isoenergetic state 
with magnetic moment of 2 [Lb only 0.003 eV higher in energy. TiNas & n d VNas clusters presented large moments 
of 3 and 4 hb respectively in their ground states. The next higher energy structures/multiplets of TiNa 5 and VNas 
are separated by significant energy differences 0.1 eV and 0.3 eV respectively from the ground states. As we discuss 
below, similar large moments are found in Ti and V doped Nag clusters and such large magnetic moments seem to 
be a feature of early TM doped small Na clusters. 

D. TMNa 6 clusters 

For TMNag clusters, we studied both planar and 3D structures having various symmetries within plane- wave USPP 
and PAW methods. In a planar geometry, we placed the TM at the center of a Na6 hexagon. In this geometry, ScNa6 
has a moment of 3 hb , while the Ti and V doped clusters have moments as large as 4 \xb and 5 [ib respectively. 

As for 3D structures for TMNag clusters, we found stable structures with C& v , and C3,, or C s symmetries as shown 
in Figure[T] In all the cases, the C5,, structure turned out to be the ground state. This pentagonal bipyramid structure 
is exactly what was found for the MgNa 6 cluster— and is also equivalent to attaching a TM atom to the Cs,, structure 
for a pure Nag cluster. 

The 3D structures turned out to be significantly more favorable in energy than the planar structure. The ground 
states of ScNa 6 , TiNa 6 and VNa 6 present moments of 1, 4 and 5 fis respectively. The next higher energy low lying 
isomers of TiNa6 and VNag have magnetic moments of 2 and 3 /ib, which are separated from the ground states by 
0.03 eV and 0.13 eV respectively. These results are presented in TableU Li et alM- have claimed that Ti and V doped 
Aug clusters have moments of 2 and 3 \xb respectively. Here we find that Ti and V doped Na clusters posses even 
larger moments. In case of V our results agree with the experimental study of Bergmann et a/ . 15 ' 16 that predict a high 
magnetic moment of 6-7 (1b for V as an impurity in bulk Na. These results, we believe, can provide an important way 
to produce magnetic clusters with large moments. In order to convince ourselves that the large magnetic moments we 
found on these clusters were not artifacts of the methods we used, and to further benchmark our calculations, we also 
calculated the properties of a CrAug cluster in a planar geometry using plane-wave USPP method for which previous 
studies exist. In agreement with Li et al^, we obtained a moment of 4 \ib- The molecular orbitals (MO) diagram 
we found was in very good agreement with that presented in Fig. 3 of ref. 13. 

One of the important questions is the nature of the polarization at the alkali metal sites. When a transition metal 
atom is embedded in a bulk free electron host, it is known that a many body state where the host polarization is 
opposed to the transition metal moment can be formed. This leads to the well known Kondo resonance^. It is then 
interesting to ask how such an effect is modified in reduced sizes. In Figure [1] we give the spin moment on the atoms 
obtained by carrying out a Mulliken population analysis of the atomic spin charges. Note that the host polarization 
changes in magnitude and sign as the cluster size is changed. These modifications can be linked to the location of the 
electronic states and their interactions between the host and the transition metal that we discuss below. 

We now try to understand the origin of such larger-than-atom magnetic moments in these clusters and the nature 
of the host polarization. Figure [5] shows the MO diagram for Na4, TMNa4 and TM atoms. The angular momentum 
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characters of the MO's are also presented. Not surprisingly, more d type MO's get occupied as the TM atom is 
changed from Sc to V. For ScNa4 in the majority a spin channel, we found that the highest occupied molecular 
orbital (HOMO), and HOMO-1 have largely d character on the Sc atom, along with some s and/or p contributions 
from the Na atoms (s,pAr a ds c in our nomenclature). In the minority f3 spin channel, the HOMO level is a mix of s 
and p type orbitals on Na (s,pjvo), and p type orbitals on Sc (ps c ). The lowest unoccupied MO's (LUMO's) in the 
a and (3 spin channels are of s,pAr a ds c type. In TiNa4 we found that the HOMO and HOMO-2 levels presented a 
s,p,/v a dn character and that HOMO-3 is mainly of dxi character. VNa4 is the first cluster to exhibit an enhanced 
magnetic moment. An analysis of the MO characters shows that the HOMO and HOMO-2 levels presented a s,pAr a 
dv character. HOMO-3 and HOMO-5 are of dy character. We can observe that in the ScNa4, TiNa4 and VNa4 series 
the spd hybridization increases, as there is an increasing number of d levels available at lower energies in the Sc, Ti 
and V series. 

Figure [3] shows the MO diagram for Na5, TMNas and TM atoms. In this series the TiNa5, VNas clusters presented 
enhanced magnetic moments. An analysis of the character of the MO's shows a spd hybridization (s,pAr a ds c ) m the 
case of ScNas. Interestingly, the MO's of TiNas and VNas present a s,pAr a Ptm (and not s,pAr a drM) character. In 
TiNas we found that the doubly degenerate HOMO, and the HOMO-1 levels are mainly of d character. This result 
shows that in addition to hybridization there is a charge transfer from the Na s levels to the Ti d states, as an isolated 
Ti atom possesses only two filled d states. The same picture was found in VNas where the degenerated HOMO and 
HOMO-2 levels are of d character. Figure Q] shows the MO diagram for Na6, TMNag and TM atoms. In this series 
also the clusters with Ti and V presented enhanced magnetic moments. Analysis of the MO's characters suggests an 
spd hybridization to be responsible for this enhancement. 

We conjecture that the major reason for the enhanced magnetic moments on these clusters is this significant spd 
hybridization. In fact, total occupancy of d and spd states in some of these clusters is larger than the occupancy 
of d states in the respective isolated TM atoms. Normally hybridization of atomic orbitals leads to quenching of 
atomic moments. Here, hybridization along with a large exchange splitting has an opposite effect. For example, 
exchange splitting (defined as the energy difference between the highest a-spin MO and the lowest /3-spin MO having 
predominantly d characters) in TiNae and VNa6 clusters is ~ 1 and ~ 1.8 eV respectively. 

A comparison of the TM doped Na6 clusters with the TM doped Au6 clusters reported by Li et al*£ is of some 
importance here. In the Aug cluster, the TM's act like guest atoms. The spectral features of the TMAu6 clusters 
showed two distinct regions which could be attributed to a TM atom and to an Aug ring. This conclusion was borne 
out by their DFT calculations which revealed that the MO's near the gap were mostly of pure e?-type on the TM 
atom, and there was a substantial gap between these levels and the deeper lying levels derived from the Au6 hexagon. 
In the present study there is no clear guest-host separation of MO's in TMNa„ clusters. On the contrary, in TiNa„ 
and VNa„ clusters, at least part of the moment comes from states that have significant spd mixing as we have argued. 

Relative values of electronegativity of Na and TM atoms may also have a crucial role to play. Na has a smaller 
(Pauling) electronegativity compared to the TM atoms. Thus, nominally one expects charge transfer from the Na to 
the TM atoms. This explains the charge transfer found in the TMNas clusters. In contrast to Na, Au has a (Pauling) 
electronegativity value much larger than those of Ti, V and Cr, suggesting that there should be a charge transfer 
from TM atoms to Au. Indeed, Torres et al. — found in their DFT calculations that the charge transfer from TM to 
Au atoms in TMAu+ clusters roughly follows the sequence of electronegativity difference between Au and TM atoms, 
and for small cluster sizes the Au-TM bond is mainly ionic. The amount of charge transfer from the TM to the Au 
atoms determines the moment on the TM atom. In TMNa„ clusters, as we have already shown, overlap of Na s and 
p and TM d type orbitals increases the moment on the TM atoms. 

The relative stability of a TMNa„ clusters can be measured by the embedding energy (EE) of the TM in the Na n 
cluster. This is defined as, 

££(TMNa„) = £(Na„) + £(TM) - £(TMNa„), (2) 

where E(TMNa„) is the ground state total energy of a TM doped Na„ cluster, E(Na„) is the total energy of a pure 
Na„ cluster, and E(TM) is the total energy of an isolated TM atom. It has been shown before that it is important 
to preserve spin multiplicity while calculating embedding energie o 35 ' 36 . Therefore, in calculating EE we used total 
energies of excited states for TM and Na n when the spin multiplicities of the ground states did not satisfy the Wigner- 
Witmer spin conservation rule 3 ^. In calculating the EE for ScNa„, TiNa4, TiNas, and VNas the ground state energies 
were employed. In the rest of the cases, excited states of either Na„ cluster or the TM atom has been considered, 
whichever gives a larger EE. The calculated EE for the TMNa n series are presented in Table [I] 

It is interesting to note certain trends in this family of clusters. For each size from Na4 to Nag, EE increases as the 
TM atom is changed from Sc to Ti and V. For any given TM atom, the EE increases from Na4 to Nas and Na6- The 
second fact suggests that in experiments, these small clusters will tend to bind with other Na atoms and form larger 
clusters. Therefore, a study extending up to larger TMNa„ clusters is of importance, and that is currently underway. 
A crucial question is, up to what size of Na clusters such large moments survive or how the moment evolves for a 
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given TM atom, as the number of Na atoms increases. This study will allow us to find out if Na clusters are able 
to present large magnetic moments in the range of 6.6 to 7 hb as found in the experimental studies of Bergmann et 
a L 15 i 16 . Our preliminary calculations show that large moments survive on the Ti and V doped Na„ clusters at least 
up to n = 922. 

In order for any cluster to be technologically useful, it has to be chemically inert. One indication of chemical 
inertness is the gap between the HOMO and LUMO levels. HOMO-LUMO gaps for these clusters are shown in 
Table U . All the clusters have a large gap of ~ 0.5 eV or more. Some of them, in fact, have gaps in excess of 1 eV. 
The large EE values and HOMO-LUMO gaps are indicative of stability in these clusters, and reinforce the idea that 
they may have useful applications. 



IV. CONCLUSION 



The present studies show that the size and structure of alkali metal clusters can significantly affect the spin magnetic 
moments of doped transition metal atoms. In particular, V doped Na4, and Ti and V doped Nas and Na e clusters 
are shown to exhibit enhanced magnetic moments. The present results for VNa„ clusters support the conclusions of 
the experimental studies by Bergman et al. 15 ' 16 that the V atoms supported on alkali layers can attain large spin 
moments. ScNag cluster also has a large moment of 3 hb in a low lying isoenergetic state. An analysis of the MO's of 
TMNa„ clusters suggests that a significant hybridization between Na s and p and TM d type orbitals coupled with a 
large exchange splitting of d type states, and a smaller value of electronegativity of Na compared to those of the TM 
atoms are responsible for the observed large moments and the host polarization. This conjecture is consistent with 
earlier calculations on neutral and cationic TMAu6 clusters. We would like to add that the host polarization can be 
expected to stabilize the orientation of the transition metal moment providing spin anisotropy with respect to the 
geometry of the alkali host. These effects are being investigated. 
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TABLE I: Electronic properties of the TMNa„ clusters. Magnetic moments (/i), HOMO-LUMO gaps (Gap), and the embedding 
energies (EE) in the ground state. Magnetic moment in the low-lying isomer, and the energy difference with respect to the 
ground state (AE). 



Cluster 




Ground state 




Low-lying 


isomer 






Gap (eV) 


EE (eV) 




AE (eV) 


ScNa 4 


1 


0.65 


1.10 


1 


0.025 


TiNa 4 


2 


1.03 


1.17 


4 


0.013 


VNa 4 


5 


0.41 


1.73 


3 


0.064 


ScNa 5 





0.51 


1.32 


2 


0.003 


TiNas 


3 


0.92 


1.50 


3 


0.104 


VNa 5 


4 


1.17 


1.89 


4 


0.283 


ScNa 6 


1 


0.67 


1.56 


1 


0.076 


TiNae 


4 


0.48 


1.89 


2 


0.030 


VNa 6 


5 


0.68 


2.21 


3 


0.128 



Figure Captions 



FIG. 1: (Color online) Ground state geometries (GS), and immediate higher energy isomers (EX) of the TMNa„ clusters. 
The bond lengths are given in Angstroms and superscripts indicate spin multiplicity. The Mulliken atomic spin charges for the 
symmetry inequivalent atoms are marked below them. 



FIG. 2: One-electron levels in Na 4 , ScNa 4 , Sc, TiNa 4 , Ti, VNa 4 and V, superscripts indicate the spin multiplicity. The 
continuous lines represent occupied levels; the dotted lines correspond to unfilled states. The degeneracy is marked by a 
number next to the level. The angular characters of the energy levels are also given. The arrows indicate the majority (up), 
and minority (down) spin states. 
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FIG. 3: One-electron levels in Nas, ScNa 5 , Sc, TiNas, Ti, VNas and V, superscripts indicate the spin multiplicity. The 
continuous lines represent occupied levels; the dotted lines correspond to unfilled states. The degeneracy is marked by a 
number next to the level. The angular character of the energy levels is also given. The arrows indicate the majority (up) , and 
minority (down) spin states. 



FIG. 4: One-electron levels in Na6, ScNa6, Sc, TiNa6, Ti, VNa6 and V superscripts indicate the spin multiplicity. The continuous 
lines represent occupied levels; the dotted lines correspond to unfilled states. The degeneracy is marked by a number next to 
the level. The angular character of the energy levels is also given. The arrows indicate the majority (up), and minority (down) 
spin states. 
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